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Status of resummed predictions for QCD final state observables 
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We provide a brief review of the current status of resummed predictions for QCD final state observables such 
as event shapes and jet rates in a variety of different hard processes. Particular emphasis is given to more recent 
developments such as the study of non global observables, development of generalized resummation formulae and 
resummations for new types of event shape variables. 



1. INTRODUCTION 

All order resummed predictions have long been 
a common method to extend the predictive power 
of theoretical QCD estimates for several observ- 
ables compared to that provided by fixed order 
pcrturbative computations alone. This applies in 
particular to observables that are measured close 
to their exclusive limits, such as event shape vari- 
ables in e + e~ annihilation 1 where the highest 
statistics are typically found close to the two jet 
limit. 

More generally, if a physical variable such 
as a typical event shape, defined for an arbi- 
trary hard process, is constrained, by observa- 
tion/measurement, to be near its Born value , 
it turns out to be sensitive to soft emissions. In 
fact if the deviation from the Born value of an 
observable (which we can take to be zero for the 
following discussion) is denoted by V, near the 
Born limit the cross-section for values up to V 
(V <§; 1) is dominated by soft and/or collinear 
logarithms. In particular the n th order perturba- 
tive estimate typically goes as : 



R™(V,Q)~ (a s (Q) In 2 V) 1 



(1) 



where the dots denote terms less singular in the 
V — > limit and Q is a scale relevant to the under- 
lying hard process. Hence the predictive power 
of fixed order perturbation theory is spoilt in the 
small V region, since the smallness of the expan- 
sion parameter a s (Q) is compensated by the pres- 
ence of large logarithms with upto two powers of 
In 1/V = L for each power of a s . 
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By now there are solid techniques in place to 
handle such logarithmic behaviour to all orders, 
which result in improved predictions valid over 
a much larger range of variable values than the 
fixed order computations alone, which become 
meaningless at very small V. These resumma- 
tion techniques rely mainly on factorisation meth- 
ods (broadly speaking) , which in turn follow from 
QCD coherence properties and the observables 
dependence on final state emissions being fac- 
torisable, usually in some appropriate transform 
space 1 . Due to the latter property not all 
variables turn out to belong to the class of re- 
summable observables-an example is jet rates de- 
fined with the JADE jet algorithm 0. 

While resummed predictions are often needed 
to improve perturbative accuracy, they also serve 
as a pathway to approach the non-perturbative 
domain one of the most challenging and certainly 
the least understood aspect of QCD. Specifically, 
to access large non-perturbative effects one has 
to enter the low k t domain, where k t is the 
transverse momentum of a typical gluon emis- 
sion, leading to a measured value of the observ- 
able. Hence in the small k t ~ VQ region, aside 
from large logarithms one also has power correc- 
tions A™/ (VQ) n where the A are non-perturbative 
coefficients. While for VQ ~ Aqcd the problem 
becomes highly non-perturbative, nevertheless for 
hard enough processes a wide range of values ex- 
ists, Aqcd <C VQ <C Q, where both perturba- 
tive resummations and non-perturbative power 
corrections are important effects. Exploring this 
region quantitatively is important as it provides 
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valuable information on the onset of confinement 
effects. 

In short resummations are a useful probe of 
QCD dynamics, which comes into its own in the 
infrared region, where parton multiplication is co- 
pious and hadronisation effects are important. In 
what follows we shall confine the discussion to ob- 
servables that do admit phase space factorisation 
and hence are resummable. We shall begin by 
considering the anatomy of a typical resummcd 
answer and indicate what is the state-of-the art 
with regard to computation of the various pieces 
we will mention. We shall then discuss some sig- 
nificant recent developments, which contributed 
to the improvement of resummed predictions and 
consequently reflect an improved understanding 
of QCD dynamics: in particular the discovery 
of non global logarithms and the advent of gen- 
eralised approaches to resummation. We shall 
also present some recent comparisons of such im- 
proved predictions with data. 

2. PROFILE OF A RESUMMED PRE- 
DICTION 

Consider an observable which has large loga- 
rithms in its perturbation expansion as below: 

R(v) = i+E < ( E R ™ ln ™ v + ° {v) ) - (2) 

n=l \m=0 ) 

One can naturally define leading logarithms 
(LL) as being the double logs m — 2n , next-to- 
leading logs (NLL) as those terms with m = 2n—l 
etc and in some cases this nomenclature is in- 
deed employed |3J- Note however that in the 
regime where a s L 2 ^S> 1, any truncation at N P LL 
order, in this terminology, is no longer useful 
since one always has neglected terms a"i m with 
In — p > m > n which are larger than one. It 
should be pointed out that the state-of-the art 
resummation enables us to go to NNLL order in 
this notation. 

However since variables that admit phase space 
factorisation typically exponentiate, it is possible 
to have a different (and in fact more commonly 
used) classification of logarithmic terms. Expo- 



nentiation means one can write [T] 

R(V) = (I + C x a s + • • -Ma s ,L) + D(V, a s ) (3) 

The first piece in brackets is simply a well be- 
haved perturbative expansion (with no depen- 
dence on V) in a s . This piece treats for instance 
the mismatch between the full real and virtual 
emission and their soft limits, included in the 
form factor E. 

All the singular dependence on V is contained 
in the form factor E. It takes the form 

VJ _ e [Lgi(a 3 L)+g 2 (a a L)+a 3 g 3 (a 3 L) + ---] 

With this exponentiation it is now possible to 
define g\ as a LL function (since it is accompanied 
by an extra logarithm L) , gi as an NLL function 
(which is purely single logarithmic a"L"), g% as 
NNLL due to the extra factor of a s relative to 
the single logarithmic g^ piece. In this notation 
the state of the art is NLL resummation in that 
resummed predictions aim to compute up to the 
(72 piece. This prediction is then matched to fixed 
order calculations. The 'remainder' piece D(V) 
vanishes in the V — > limit . 

Having clarified the basic notation and termi- 
nology we now turn our attention to some more 
recent developments. We begin the next section 
by providing a short discussion on non global ob- 
servables. 

3. NON GLOBAL EFFECTS 

Until very recently the resummed predictions 
available in the literature were all made employ- 
ing an independent emission formalism. In other 
words the matrix element for multiple gluon pro- 
duction was approximated by a factorised prod- 
uct of single soft emission contributions: 

dV n « Y\ C F a s (ku)dr/— —. (5) 

1=1 

with k t the transverse momentum and rj the ra- 
pidity wrt the emitting jet axis. 

The independent emission pattern (in one form 
or another) was used in conjunction with phase 
space factorisation to yield the form factor E. In- 
dependent emission essentially means one can as- 
sume each final state gluon to be emitted directly 
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by the hard initiating parton as in an abelian the- 
ory. Subsequent branching of these gluons can in 
fact be neglected for several observables, apart 
from the contribution of such decays to the run- 
ning of the coupling. 

However the independent emission approxima- 
tion © was found to be insufficient at single log- 
arithmic level (ie at the level of the function c/2) 
for several observables 0]. In some of these cases, 
such as energy flow away from jets single loga- 
rithms were in fact a leading effect due to absence 
of collinear enhancements, which meant that the 
function g\ was absent. In general the observables 
for which independent emission breaks down at 
single-log level have one common feature respon- 
sible for this breakdown - they are sensitive to 
emissions only in a limited rapidity region. For 
this reason such observables and corresponding 
logarithms are refered to as non global, while the 
complementary set of observables is refered to as 
global. 

The origin of the problem is simple to under- 
stand and briefly explained below. 




picted above. Then by restricting, for example, 
the transverse energy in Q such that Qq <C Q 
with Q a hard scale, one is vetoing real radiation 
above scale Qq into Q. Assuming an indepen- 
dent emission ansatz one needs to just veto direct 
emission from the hard jet lines as in diagram (a) 
above. 



10 L 



10" 



10" 



10" 



10" 







V \ 
V \ 
•N \ 








(a) 




w 

V 

\ 


\ \ 












^ \ 
\ \ 

\ \ 

\ 

\ 

\ 






- Pate 


;h with Ar| = Ac 
3 with Ar| = 1 
ite slice 


) = 2 


\ \ 
*\ \ 
^ \ 
\ \ 




-- Slic 
- Infir 




\ 1 

\ 



10" 

0.1 0.2 0.3 0.4 0.5 0.6 
t 

Figure 1. Non global function S(t) plotted 
against t = Jq^ ^-a s (k t ) for different definitions 
ofO 



This yields a single logarithmic form factor 
E SL ~ exp[-a s In Q/Qq] . (6) 



Consider a measurement of say energy flow 
away from hard jets in a solid angle O as de- 



The above result is 
This can be seen by 
comparing with the 
an exact fixed order 
an additional source 
starting from 0{a 2 s ) 
ancy with the single 
estimate. 



incorrect at single log level, 
expanding it to 0(a 2 s ) and 
logarithmic dependence of 
estimate at 0{a 2 s ). In fact 
of single logarithms is found 
which removes the discrep- 
logarithms in the fixed order 
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At order 0(a 2 s ) this is just emission of a soft 
gluon with energy u> ~ Qq by a relatively harder 
gluon with energy to' with Qq <C <C Q outside 

(see diagram (b)). This and analogous higher 
order contributions are missed by the indepen- 
dent emission approach which would be correct if 
both the gluon outside fl and the softest gluon in- 
side f2 were measured (vetoed) . The effect of the 
softest gluon would cancel against virtual correc- 
tions at single logarithmic level. However since 
we are only measuring inside SI the softest gluon 
emission is significant even in the presence of the 
harder emitter outside f2 and hence the gluon 
branching contribution appears. 

The generalization of the above effect to n th 
order is the coherent emission of a single soft- 
est gluon into f2 by an ensemble of n — 1 glu- 
ons outside fl (diagram (c) above. The n — 1 
gluon emitters are themselves ordered in energy, 
loq u)i oj n -i . The effect of such multiple 
wide-angle (non collinear-enhanced) soft emission 
is precisely a single logarithm, a"ln™^-. Re- 
summing these terms needs a change of approach 
from the independent emission approximation. In 
Ref. ' 115) , the resummation of single logarithms is 
carried out numerically by using the dipole evo- 
lution picture which captures the essential sin- 
gle logarithms in the large N c limit. The result- 
ing function is plotted in Fig. ^ Following this 
a non-linear evolution equation corresponding to 
the dipole evolution was derived in Ref. [S] and 
its numerical solution yielded identical results to 
those obtained in Ref. 

Subsequent to the discovery of non global logs, 
there was some effort made at understanding how 
to define observables in a way such that non 
global effects may be minimised, which would 
reduce any uncertainty to do with missing non 
global 1/N 2 suppressed effects which cannot yet 
be resummed. The reduction or removal of the 
non global component is also useful for facilitat- 
ing the phenomenology of multi jet event shapes 
and energy flows. This is because non global logs 
have thus far only been explicitly computed in 
the two jet case, although a similar structure is 
expected in the extension to multi-jet events and 
their calculation will also be along the same lines, 
employing dipole (non-linear) evolution. 



In the above regard Appleby and Seymour con- 
sidered the effect of jet clustering algorithms on 
non global logarithms. They found that if they 
defined rapidity gaps in terms of minijet en- 
ergy flows (with the minijets being defined by a 
clustering alogorithm), rather than a sum over 
hadronic energies in the gap, the non global com- 
ponent was significantly reduced. This is essen- 
tially because the clustering algorithm has the ef- 
fect of pulling soft hadrons (partons) out of the 
gap and clustering them with harder emissions 
outside the gap. Hence the non global contribu- 
tion is not triggered except in specific geometrical 
configurations which survive the clustering, which 
reduces its numerical significance (see \7\ for de- 
tails). 

Berger, Kucs and Sterman introduced event- 
shape/energy flow correlations [Hj aimed at con- 
trolling non global effects in energy flow. This 
meant simultaneously restricting the energy Qq 
in fl alongside limiting the value of an event shape 
variable V, defined in terms of hadron momenta 
outside O. Doing so amounts to restricting trans- 
verse momenta of soft gluons both outside the gap 
region, by controlling the event shape k t < VQ, 
and inside the gap k t < Qq with Qq < VQ <§C Q. 
Non global or secondary logarithms now appear 
as a™ In" contributions and hence varying the 
values of V and Qq in tandem, allows one to con- 
trol the significance the of the non-global contri- 
bution. 

In the course of these studies the same authors 
also introduced new kinds of event shapes with an 
adjustable parameter a that allows one to control 
the approach to the two jet limit [HJ. These vari- 
ables have a parametric dependence on the trans- 
verse momenta k t and rapidity r\ with respect to 
the thrust axis, of the form V = J2t kue~ {1 ~ a)lv ' 1 , 
where the sum extends over all final state partons. 
The resummation to NLL accutacy, for these ob- 
servables can be found in Ref. [S]. A special case, 
a = 0, of the above class of variables is the much 
studied thrust variable. 

Dokshitzer and Marchesini [H] further extended 
the study of event shape/energy flow correlations 
by showing that the non global part of the an- 
swer, which involves logarithms of QV/Qq, fac- 
torises from the usual global resummation of the 
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event shape V (this time defined as a sum over 
all final state hadron contributions in and outside 
the gap). 

Although some effort has been devoted to 
understanding and computing non global loga- 
rithms, the only inclusion of this effect in com- 
parisons to data thus far has been in the case of 
DIS event shape variables. For details the reader 
is refered to Ref. 10 . A recent comparison of 
a matched resummed prediction (including non 
global effects and power corrections) to HI data 
is shown in Fig. 2, for the thrust defined wrt 
the thrust axis in the Breit current hemisphere 
of DIS. 



4. THREE JET SHAPE VARIABLES 

An interesting and fairly recent development 
is that of NLL resummed predictions for three 
jet event shapes in e + e~ , DIS (two final state 
+one incoming jet) and hadron-hadron collisions 
(two incoming jets and a final state (radiated) 
jet+vector boson). For an example see (TT]. For 
all these processes, at Born level one finds that 
the three hard jets lie in a plane defined say by the 
thrust and thrust major axes. Small deviations 
from the Born level can be accesed by studying 
small values of the out of plane momentum K out . 
Compared to two jet event shapes such as thrust, 
the three jet variables offer some new challenges. 
From the point of view of the resummed predic- 
tion, one as to also take account of coherent, soft 
interjet radiation while in the global two jet case 
only radiation emitted collinearly (effectively in- 
coherently) by each jet contributes to NLL accu- 
racy. The interjet radiation gives a geometry de- 
pendence at NLL level, characteristic of three jet 
variables. Such variables would be very interest- 
ing to study experimentally especially from the 
viewpoint of hadronisation corrections. This is 
a particularly important intermediate step along 
the way to understanding power behaved hadro- 
nisation effects in multijet configurations, such as 
in dijet event shapes in hadroproduction (2 — > 2 
processes) . 
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Figure 2. Resummed distribution for the thrust 
defined wrt the thrust axis in the current hemi- 
sphere of the DIS Breit frame compared to HI 
data. From top to bottom represents increasing 
Q values from 15 - 81 GeV. 



5. GENERALISED FORMULAE AND 
AUTOMATED RESUMMATIONS 

Of late there has been significant effort and 
progress made in developing resummation formu- 
lae that are as general as currently possible and 
which reduce the need for treating different ob- 
servables and different processeses on a laborious 
case-by-case basis. The idea is to write down 
master formulae which have general applicability 
and contain process and observable dependence 
in a few parameters which need to be computed 
on a case specific basis [12113114] , For example 
processes involving an arbitrary number of hard 
partons were treated in Ref. 13 . 

A generalised approach that aims at facilitat- 
ing the computablity of resummed results , us- 
ing numerical methods, and in fact automating 
the entire process of resummation was developed 
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in The idea roughly, is to test the observ- 

able's dependence on soft and collinear emissions 
and determine numerically the parametric depen- 
dence on transverse momenta, rapidity and az- 
imuth wrt the nearest leg. Once this depen- 
dence is obtained, the results are automatically 
inserted into a master formula valid for an arbi- 
trary case (there however being some conditions 
such as globalness, that the observable has to sat- 
isfy for the master formula to be valid). Final 
results are numerically obtained by means of a 
computer program. This is a particularly pow- 
erful method since almost no analytical effort is 
needed to generate the final results and the entire 
procedure is automated. Its advantage is partic- 
ularly manifest in cases where analytical calcu- 
lations are too cumbersome, requiring the inver- 
sion of multiple Mellin and/or Fourier transforms 
or the observable depends in a complicated way 
on multiple soft emissions. A good example of 
this latter point is the Durham jet finding algo- 
rithm for which the analytic computation of the 
full NLL function g-i was intractable due to effects 
to do with recombination of soft partons. The nu- 
merical resummation methods were on the other 
hand able to determine this function for the first 
time . 

Another example of the power of the auto- 
mated resummation approach was provided by 
the computation, for the first time, to NLL ac- 
curacy of an event shape in hadron-hadron col- 
lisions (with a dijet final state) - the transverse 
thrust distribution This is defined as follows: 



W±i-n±\ 



(7) 



where as usual the sum extends over all final state 
particles and p± is the transverse momentum wrt 
the beam direction while ft is a unit vector found 
by maximising the sum in the numerator. The 
results for the different dijet production channels 
is illustrated in Fig. 2. 

In conclusion it is worthwhile to note that 
there has been significant progress in the develop- 
ment of resummed computations in the last few 
years. The consequences of these developments 
are naturally of great value to experimental stud- 
ies and future QCD phenomenology. Addition- 



ally they point towards unravelling of previously 
unexplored QCD dynamics and therefore to im- 
portant progress on the theoretical front. 
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Figure 3. Plot of the transverse thrust distribu- 
tion in hadronic collisions at y/s = 1.96 GeV.. 
Fig. taken from Ref. |14| . 
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